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Recently it was argued that it might be possible to treat the conventional nuclear 
structure problem - nonrelativistic point nucleons interacting through a static and rather 
singular potential - as an effective theory in a shell-model basis. In the first half of 
this talk we describe how such a program can be carried out for the simplest nuclei, the 
deuteron and 3 He, exploiting a new numerical technique for solving the self-consistent 
Bloch-Horowitz equation. Some of the properties of proper effective theories are thus il- 
lustrated and contrasted with the shell model. In the second half of the talk we use these 
examples to return to a problem that frustrated the field three decades ago, the possibility 
of reducing the effective interactions problem to perturbation theory. We show, by exploit- 
ing the Talmi integral expansion, that hard-core potentials can be systematically softened 
by the introduction of a series of contact operators familiar from effective field theory. The 
coefficients of these operators can be run analytically by a renormalization group method 
in a scheme-independent way, with the introduction of suitable counterterms. Once these 
coefficients are run to the shell model scale, we show that the renormalized coefficients 
contain all of the information needed to evaluate perturbative insertions of the remaining 
soft potential. The resulting perturbative expansion is shown to converge in lowest order 
for the simplest nucleus, the deuteron. 



1. INTRODUCTION 

It is an honor to be able to address a distinguished audience at this celebration for 
Achim Richter, a very good friend and valued colleague who has done so so much to 
advance nuclear physics. It is also great fun to return, in this talk, to a topic that amused 
several of the theoreticians in this room in the 1970's, the possibility of a perturbative 
expansion of effective interactions and operators. The thesis of this talk is that a mar- 
riage of the formidable technology of the shell model (SM) with the modern ideas of the 
renormalization group and effective theory (ET) may help us make some progress on this 
important problem. 

In the first half of this talk we contrast the conventional SM with the results of a proper 
ET carried out in a SM basis. We introduce a new and rather efficient numerical method 
for solving the self-consistent Bloch-Horowitz equation governing effective interactions 
and operators. The method, while numerical, systematically extracts from the excluded 
high-momentum space the information needed to construct effective interactions, and 
thus is in the spirit of ET techniques. The results for two simple nuclei, the deuteron and 
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3 He, illustrate many of the properties of correct effective theories (simple wave function 
evolution in the included space, nontrivial normalizations) that are absent in the shell 
model. 

The effective interactions problem is known to be highly nonperturbative. In the second 
half of this talk we exploit these solutions in an effort to understand the source of the 
nonperturbative scattering at high momenta. The first source we encounter arises from 
the overbinding of harmonic oscillator wave functions. A straightforward reorganization 
of the Bio ch- Horowitz equation - equivalent to resumming the kinetic energy operator 
to all orders - removes this problem. The next source is hard-core scattering. To attack 
this problem we build a bridge between effective field theories and potential problems 
by showing that the contact operators introduced in the former generate a version of 
the Talmi integral expansion in a SM basis. Thus a rather singular NN potential can be 
systematically softened by removing the lowest order Talmi integrals, replacing these with 
contact potentials. The effective interactions problem for these contact interactions can 
be solved analytically via a shell-by-shell renormalization group equation. The hard-core 
effective interactions problem corresponds to understanding the running of the coefficients 
of these contact interactions. The running can be done in a scheme-independent way - that 
is, the effects of the hard core on every matrix element in the SM space can be evaluated 
exactly - with the introduction of a finite number of higher-order contact interactions (or 
counterterms) . 

With the hard core problem thus solved, we return to the remaining soft parts of 
the potential, showing they can be inserted perturbatively among the infinite hard-core 
scattering series. The only information needed is the renormalized coefficients of the 
contact interactions. We illustrate, for the lowest order calculation in the deuteron, that 
this procedure reduces the effective interactions problem to perturbation theory. 

Although the talk focuses on the simplest possible example, the techniques employed 
are quite general. Details not presented here can be found in [0-0]. 

2. EFFECTIVE INTERACTIONS 

Consider a collection of nonrelativistic nucleons interacting through a potential, e.g., of 
the Nijmegen or Argonne/Urbana type, within an infinite Hilbert space 

H = \ £(Ty + V«,-), (1) 

where T^- is the relative nonrelativistic kinetic energy operator and V%j the nucleon-nucleon 
potential. The related SM Hamiltonian acts in a restricted space and employs a softer 
"effective" potential, 

HsM = \j:{T l3 + Vf). (2) 

Motivating Hsm is the notion that the determination of V e ^ might be simpler than 
solving the original A-body problem: the foundation of Brueckner theory is that high- 
momentum contributions to the wave function might be integrated out in a rapidly con- 
verging series in p nuc iear or, equivalently, in the number of nucleons in high-momentum 
states interacting at one time outside the SM space. 
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As will become apparent later, a typical SM space will contain explicitly ~ 60% of the 
wave function that resides at long-wavelengths, thus guaranteeing at least a qualitative 
description of A-body correlations important to soft collective modes. The strength of the 
SM is its capacity to treat this long-wavelength part of the problem nonperturbatively: 
the technology for direct diagonalizations in large SM spaces is quite remarkable, including 
recent progress in Lanczos-based methods M, in treatments of light nuclei involving many 
shells 0, and in Monte Carlo sampling |||7]|. Unfortunately the SM - it is a model - 
makes a number of uncontrolled approximations that distinguish it from an ET. These 
are discussed in Ref. jl] and include effective interactions lacking the correct functional 
form; a misunderstanding of the proper normalization for SM wave functions; dependence 
of results on fictitious parameters such as starting energies; and the almost total neglect 
of effective operators. 

2.1. Self-consistent B loch- Horowitz Solutions 

The effort to reformulate the SM as an ET begins with a definition of the "SM" space. 
The goals of handling bound states and of generating an effective interaction that is trans- 
lationally invariant leaves one sensible choice: many-body states constructed from har- 
monic oscillator Slater determinants. To exploit the relative/center-of-mass separability of 
harmonic oscillator Slater determinants, one must separate the SM and high-momentum 
spaces so that all configurations satisfying 

E < Asm^sm (3) 

are retained in the former. For example, a SM calculation of 16 with Asm — 4+A , where 
A is the number of quanta in the 16 O closed shell, would include all ATiujsm configurations, 
e.g., OpOh, 2p2h, and 4p4h excitations of nucleons from the lp shell into the 2sld shell, 
Iplh excitations of a Is shell nucleon into the 3s2dlg shell, etc. One can define the 
projection operator onto the excluded high-momentum space by 

QSM = Q{A S M,b S A4). (4) 

where bsM is the oscillator parameter - defined here in the usual SM way for independent 
particle motion - and Uujsm = tjti — is the corresponding energy. Thus the included or 

1V10 SM 

"SM" space is defined by two parameters, Asm an d bsM ■ The preservation of translational 
invariance is also important numerically, as it reduces the two-body ladder to an effective 
one-body problem, etc. 

The resulting Bloch-Horowitz (BH) equation |§ is then 

H*" = H + H 1 Q SM H 

H eff \V S M) = E\S>sm) \^sm) = {1-QsmW) (5) 

where is the exact wave function and H\fy) = E\^). The difficulty posed by this 
equation is the appearance of the unknown energy eigenvalue in the equation for H e ^ . 
Thus this system must be solved self-consistently. Note that there is no explicit reference 
to the harmonic oscillator in this equation: it enters only implicitly through Qsm i n 
distinguishing the long-wavelength "SM" space from the remainder of the Hilbert space. 



4 



We emphasize that a proper solution of the BH equation must yield results (energies and 
operator matrix elements) that are independent of Qsm- 

A new procedure for solving this problem was introduced in Ref. jl]. It is based on a 
mapping of the full Hamiltonian in the high momentum space - this space can be made 
finite, extending to some scale Aoo ~ 3 GeV adequate to fully resolve the hard cores of 
realistic potentials - into a simpler, truncated Hamiltonian via the Lanczos algorithm, 
followed by the demonstration that the the BH Green's function can be constructed from 
the Lanczos matrix as a function of E, thus solving the effective interactions problem. 
That is, one can iterate on E simply by repeating a shell model calculation, as the effective 
interaction can be immediately generated for any desired E. In practice the procedure 
converges very rapidly. 

In our view the Lanczos approach to the effective interactions problem appears to be 
remarkably simpler than the standard procedures of the field, particular in view of the 
need for extensions to multi-nucleon ladders. The traditional approach divides the ef- 
fective interactions problem into an energy-independent piece (often called the Q-box) 
and an energy-dependent one (represented by folded diagrams). Frequently the energy 
dependence is removed by making a unitary transformation to a non-Hermitian effective 
Hamiltonian. The Lanczos procedure is also much more in the spirit of EFT: the full 
Hamiltonian in the high-momentum space is never constructed. Rather, it is replaced by 
a much smaller Lanczos matrix which contains exactly the most relevant long- wavelength 
information of the full matrix, the 2n-l lowest moments (here n represents the nth itera- 
tion in the Lanczos algorithm). Thus the procedure can be viewed as a numerical ET in 
which this information is recursively extracted. Self-consistent solutions must be obtained 
for each state, thus yielding a Hermitian but energy- dependent effective interaction. 

2.2. Properties of Effective Wave Functions and Operators 

A first test of the techniques outlined above is to solve the BH equation for some SM- 
like space to then see if the resulting self-consistent energy is, indeed, the correct value 
(e.g., in agreement with Faddeev or other exact calculations). For model spaces of 2, 4, 
6, and 8Hujsm in the case of the deuteron we obtained a binding energy of -2.224 MeV 
(using V2b SM = 1.6/ and A^ = 140). The exact result is -2.2246 MeV. Our 3 He result 
for same SM spaces and A ro = 60 is -6.87 MeV, in agreement with the corresponding 
Green's function Monte Carlo (GFMC) result of -6.87 ± 0.03 MeV. Note that energies 
are variational in A^, decreasing as A^ is increased. 

More interesting is the evolution of the wave functions, shown in Table [I], which is quite 
unlike that of typical shell model calculations. The wave functions obtained in different 
model spaces agree over overlapping parts of their Hilbert spaces. Thus as one proceeds 
through 2hu>sA4, 4huJsM, G^sm, ••• calculations, the ET wave function evolves only by 
adding new components in the expanded space. The normalization of the wave function 
grows accordingly. (This normalization must be calculated, as described below.) Thus, 
for 3 He, the OHcusm ET calculation contains 0.311 of the full wave function in the effective 
space; the 0+2+4husM result is 0.700. 

Although there is clearly an intimate relation between effective interactions and effec- 
tive operators, it is standard in the shell model to calculate nuclear responses with bare 
operators, or perhaps with bare operators renormalized according to effective charges de- 
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Table 1 

ET results for the 3 He ground state wave function calculated with the Argonne vlS po- 
tential. The columns on the right correspond to different choices of the ET model space, 
the analog of a SM space. The rows correspond to the resulting amplitudes for the des- 
ignated, selected configurations. The quantities within the parentheses are the square of 
the norm of the effective wave function, e.g., the fraction of the 3 He ground state that 
resides in the "SM" space. The basis states are designated somewhat schematically as 
| N,a), where N is the total number of oscillator quanta and a is an index representing 
all other quantum numbers. 



state 


amplitude 


OkuJsM 


2%UJSM 




QflUJsM 


8 husM 


exact 


(31.1%) 


(57.4%) 


(70.0%) 


(79.8%) 


(85.5%) 


(100%) 




0,1) 


0.55791 


0.55791 


0.55791 


0.55795 


0.55791 


0.55793 




2,1) 


0.00000 


0.04631 


0.04613 


0.04618 


0.04622 


0.04631 




2,2) 


0.00000 


-0.48255 


-0.48237 


-0.48243 


-0.48243 


-0.48257 




2,3) 


0.00000 


0.00729 


0.00731 


0.00730 


0.00729 


0.00729 




2,4) 


0.00000 


0.16707 


0.16698 


0.16706 


0.16706 


0.16708 




4,1) 


0.00000 


0.00000 


-0.02040 


-0.02042 


-0.02043 


-0.02047 




4,2) 


0.00000 


0.00000 


0.11267 


0.11274 


0.11275 


0.11289 




4,3) 


0.00000 


0.00000 


-0.04191 


-0.04199 


-0.04208 


-0.04228 




4,4) 


0.00000 


0.00000 


0.28967 


0.28978 


0.28978 


0.29001 



termined phenomenologically at q 2 = 0, using SM wave functions normed to 1. As we 
now have a series of exact effective interactions corresponding to different model spaces, 
we can test the validity of this approach. The results for the elastic magnetic form factors 
for the deuteron and 3 He are shown in Fig. 1. 

One sees in each case that by the time one reaches a momentum transfer q ~ 2.5//, 
random numbers are being generated: bare operators used in conjunction with exact 
effective wave functions generate results that differ by an order of magnitude, depending 
on the choice of the model space. This is not surprising: if one probes the nucleus at 
momentum transfers > 2kp, where hp is the Fermi momentum, most of the resulting 
amplitude should reside outside the long-wavelength model space. That is, the strength 
resides entirely in the effective contributions to the operator. If these components are 
ignored, the results have to be in error. 

Clearly the effective interaction and effective operator have to be treated consistently 
and on the same footing. If O is the bare operator, one finds 

(*/|0|*i> = {^ e / f \O eff \^1 ff ) (6) 

where 

= (1 + HQ SM - L— - )6(1 + - 1 Q SM H) (7) 

Hi f — II LJsM hji — LisMli 
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and where the effective wave function normalization of \^ ) and ), mentioned 
earlier, must be determined using the effective operator 1, e.g., 

1 = = (^ e// l(l + HQsm w l n )(1 + p ^ ttQsmH) I ^ e// ) (8) 

These expressions can be evaluated with the Lanczos Green's function methods described 
earlier. When this is done, all of the effective calculations, regardless of the choice of the 
model space, yield the same result, given by the solid lines in Fig. 1. 

We would argue, based on this example, that many persistent problems in nuclear 
physics — ranging from the renormalization of Qa in (3 decay to the systematic differences 
between measured and calculated Ml electromagnetic form factors — very likely are due 
to naive treatments of operators, treatments that fail to satisfy the basic rules of ETs. It 
should be apparent from the above example that no amount of work on H e ^ will help 
with this problem. What is necessary is a diagrammatic basis for generating H e ^ that 
can be applied in exactly the same way to evolving O e ^ . From this perspective, phe- 
nomenological derivations of H e ^ by fitting binding energies and other static properties of 
nuclei are not terribly helpful, unless one intends to simultaneously find phenomenological 
renormalizations for each desired operator in each q 2 range of interest. 

The point of this discussion has been largely pedagogical. To our knowledge the exam- 
ples given above are the only ones in the classical nuclear physics literature in which a 
"model space" calculation has been formulated in a way that satisfies the basic rules of 
ET. The results — particularly the invariance of energies and operator matrix elements 
under changes in how we define Q — are obvious from the perspective of ET. The fact 
that standard techniques like the SM violate so many of the rules of ETs is cause for 
optimism: much can be done to improve the rigor of such nuclear physics tools. As we 
have demonstrated here, new numerical methods can be developed to handle the full ET 
problem at the cost of only modest additional effort. 

3. RESUMMATIONS: LONG AND SHORT DISTANCES 

The effort to derive SM effective interactions directly from realistic NN interactions was 
largely abandoned in the early 1970's when it became apparent that the high-momentum 
scattering problem was nonperturbative. A series of efforts involving perturbative expan- 
sions in either the potential V or the corresponding G-matrix (the sum of the two-body 
ladder diagrams) yielded discouraging results, with higher-order terms often dominat- 
ing over lower-order onesflO]. In particular Shucan and Weidenmuller showed that any 



overlap in the spectra of included and excluded, high-momentum states would lead to a 



nonperturbative series [11,12] . Thus, the uncontrolled approximations in the SM approach 
grew out of necessity: it became a largely phenomenological tool in part because early 
efforts to provide a more rigorous basis for the model ended in failure. 

Given our two SM-like ET calculations, the deuteron and 3 He, it is interesting to 
return to this old problem of the nonperturbative behavior of the nonrelativistic many- 
body problem, using these examples as laboratories. Today's challenge, in fact, is more 
modest than that of the 1970's: modern computers coupled with new algorithms, like the 
Lanczos solution of the BH equation introduced above, allow us to solve rather complex 
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Figure 1. The magnetic elastic form factor for the deuteron (top) and 3 He (bottom) cal- 
culated with the exact H e *' , SM wave functions normalized to unity, and a bare operator 
are compared to the exact result (solid line) . When effective operators and the proper 
wave function normalizations are used, all results become identical to the solid line. 



nonperturbative problems. To extend the calculations above to 4 He and heavier nuclei 
one has to achieve only modest increases in efficiency, such as a substantial lowering of 
Aqo. Below we summarize some of our most recent work || that suggests that the goals 
of the 1970's — a perturbative treatment of the nonrelativistic nuclear physics problem 
- may not be out of the question. Through renormalization group and other techniques, 
the most severely nonperturbative contributions to the potential can be systematically 
extracted and treated. In effect, we try to show how many of the ideas of EFT can be 
adapted to the nuclear potential problem, without abandoning the basic idea of the SM 
that rigorous nonperturbative diagonalizations in a long-wavelength included space are 
important to the physics. 



8 



3.1. Harmonic Oscillator Overbinding and Kinetic Energy Resummations 

We begin by studying the behavior of SM effective interaction matrix elements (a\H e f* \/3 
in a perturbative expansion of the BH propagator 

1 1 + (9) 



E — QH E — Hq E — Hq E — Hq 

evaluating the series term by term. Here H Q and Vq are the harmonic oscillator Hamilto- 
nian and potential. We thus lower Aoo and test whether states above A^ can be treated 
perturbatively. 

Once a scale of ~ 70 is reached, striking differences appear in the rate of convergence 
of different matrix elements. All are greatly improved in first and second order perturba- 
tion theory. Those matrix elements well below the "boundary" at A OT converge quickly 
to their correct values in this way. However, those matrix elements (a\H e ^\/3) where 
either |a) or \(3) resides in the last shell (with N = A^) converge only very slowly to the 
correct values, after the initial improvement in low-order perturbation theory. Typically 
~ 10 3 orders of perturbation theory are required to produce the correct value. This clearly 
suggests that slow convergence is associated with the relative kinetic energy operator QT 
contribution to QH, as the only transitions to states outside the Hilbert space generated 
by this operator have An = 1. At large r the strength of this transition becomes quite 
large, —QVo(r)\a), reflecting the unphysical asymptotic behavior of harmonic oscillator 
wave functions. This amplitude propagates nonperturbatively as (V — Vo)/H ~ 1. Even- 
tually enough high-momentum harmonic oscillator wave functions are coupled together 
to produce the softer asymptotic fall-off characteristic of the correct bound-state wave 
function. 

As discussed in Ref. [pj, the required resummation is guided by the observation that 
the true potential, V(r), falls off properly at large r. Thus a reorganization of the BH 
equation in which the propagator is always sandwiched between V(r) should remove the 
unwanted propagation. This leads to the following recasting of the BH equation: 

(a\H^\P) = (a\T\(3) + ((a\-(a\)E-T(0)-\(3)) 

+(a\V + V E ^ QH QV0) (10) 

where 

la) = ^|a> (11) 

If a and (3 are not in the last shell, \&) = \a) and \(3) = so that the above rewriting 
of the BH equation reduces to the original form. Otherwise, a modified wave function 
generated by Eq. ( |TT|) must be used. This rather simple resummation removes the worst - 
though not the most interesting - source of nonperturbative behavior in the BH equation. 

It should be clear that a large reduction in the scale A^ at the cost of an evaluation of 
Eq. ( |TTD represents a tremendous savings: in place of a dense matrix QH whose elements 
have to be evaluated numerically, we have a sparse matrix QT whose matrix elements are 
known analytically. For example, in the case of the deuteron, QT is triadiagonal in the 
harmonic oscillator basis, allowing us to write the Green's function in a continued fraction 
expansion. 
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3.2. Talmi Integrals, Renormalization Group, and the Hard Core 

Once the replacement |a) — > \a) is made, can be lowered to ~ 40 while maintaining 
1 keV accuracy, if one works to third order in perturbation theory. But with further 
lowering, errors arise in the perturbative expansion that persist even if calculations are 
carried to very high order. At Aoo = 30, third order perturbation theory reduces ~ 
10% errors in (n' = IV = 0\H e ^\n = 11 = 0) to ~ 0.2%. But an error in excess of 
0.1% — corresponding to 50 keV — persists after 10 additional orders of perturbation. 
Numerically one can verify that the nonperturbative tail is generated by the scattering 
at very small r. As expected, the nonperturbative contributions to s — d matrix elements 
are much smaller than those for (n' = IV = 0\H e ^\n = 1/ = 0) and other s — s matrix 
elements. 

We think this hard-core problem can be solved. To keep the presentation reasonable, 
below we concentrate primarily on the lowest order (LO) treatment. (The NLO and 
NNLO results are given in Ref. 0.) One begins with the Talmi integral expansion for 
radial harmonic oscillator matrix elements in the relative coordinate 



(n'l'\V sr (r)\nl) = / R nll i(r)V sr (r)R nl ( r ydr 
Jo 

= ^£«Z';n,Z;p)/ p (&) (12) 



where 



Kb) - 2 rV 2/fe V (r) r2P+2dr (13) 

Ip[b) ~ T(p + 3/2) Jo sr{ ) ' [L6) 



b = V2b 



SM- 



The B(n',l';n,l;p) are known analytically ||13||. Here V sr (r) denotes the short-range con- 



tributions to the potential (we return to this point below). The integrals I p can be viewed 
as a systematic expansion of the nonperturbative hard core in terms of the parameter 
(r c /b) 2 , where r c is a distance associated with the size of the hard core that remains 
unresolved at scale A^. Thus the LO term in the expansion is 

*=wf w * (14) 

while in NLO one includes 



42. = rm&L VJrydr. (15) 
In NNLO one obtains 



^ = ww^C Wdr - (16) 
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and so on. 

Now the radial component of the s-wave potential Vo<5(r) is 

Vf\r) = v ^- 2 5{r). (17) 



Substituting this into Eq. fll2"D yields 

(n'l' = 0\V 5 (0) (r)\nl = 0) = B{n\ 0; n, 0; 0)Iq(&) 



where 



™ = rMirf"' < 19 > 

It follows that the contact interaction will produce the exact lowest-order Talmi integral 
contribution (Eq. (|TiD ) to V sr provided 



oo 

2, 



V = 4ir V sr (r)r 2 dr. (20) 
J o 

Using a normalization that will be convenient later, we can then rewrite V sr (r) as 
Vsr(r) = V<f>(r) + 6 3 y^a-(A oc )5(r) (21) 

where u = ^ and Vjfi is a new potential whose leading-order behavior is of order (r c /b) 2 
relative to V sr (determined by subtracting from the original Talmi integral expression for 
V sr the contribution from Eq. flHD). The coefficient a^A^) is a dimensionless coupling. 
(Note that throughout this discussion, we have suppressed the spin and isospin of Vg r (r). 
Thus Vsrir) represents the radial function obtained after spin and isospin matrix elements 
have been taken. For example, the components of the avl8 potential contributing to s — s 
transitions is 

V^t) = V^r) + V 2 ir)n ■ r 2 + V z (r)B x ■ a 2 + V A {r)a t ■ a 2 r x • r 2 (22) 
so that 

V S rir) = Vir) - 3V 2 (r) + V 3 (r) - 3V 4 (r) (23) 

for the deuteron (1 — 0, s — 1, t — 0).) 

The above procedure is, in fact, general. The NLO contribution (Eq. (|15|)) can be 
removed from s — s wave deuteron matrix elements of Kr( r ) by introducing a second, 
higher-order contact operator 



s (Aoo)i(V 2 5(r) + ,5(r)^ 2 ) (24) 

where 



a 2 [--w/ 2 ' 



a?(Aoo) = L VUrydr, (25) 
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while the NNLO contributions can be removed by 

3 

10 



af (Aoo)(V 2 5(r)^ 2 + + 5(r)^ 4 )) (26) 



where 

<(Aoc) = / V sr (r)r e dr. (27) 

yno'Tiuj Jo 

Likewise there are coeficients of additional contact operators (denoted a s 2 d , a| d , and af d in 
Ref. ||) that remove the Talmi integral contributions to s — d and d — d matrix elements 
to NNLO. 

Up to this point we have simply rewritten the "bare" potential — the potential that 
acts in the space defined by Aoo — in an entirely equivalent form, exploiting the Talmi 
integral expansion and the fact that the unique translationally-invariant two-body contact 
operators generate those Talmi integrals. Now we consider, in LO, the effects on a^ 3 of 
integrating out the single shell at A^, thereby mapping the original problem into an 
effective one with A^ — > A ro — 2. The resulting BH equation instructs one to simply 
renormalize the coefficient a^ s (A) in LO H, 



r(A±3) a s S{A y 



(A - 2) = <(A) + (2f 



2)- E -(A + l)-a s °(A) 



2/ ~U V-y (A)] 

where E = E/hu is a. dimensionless energy. (£" is the BH energy being determined self- 
consistently. ) This leading-order RG equation is an operator equation, properly correcting 
every matrix element in the new space below A^ — 2. Clearly this difference equation can 
be used to run ao s (A oc ) to any desired new scale, including Asm- 

While the RG equation is a difference equation, the connections with EFT become even 
clearer by introducing the natural definition of the derivative 

<*ag' (A) _ aS'(A)-qg'(A-2) 

-^A~ = 2 • (29) 

As will be apparent from results given below, the terms in the denominator of Eq. ( p8|) 
that compete with A are one power lower in A. Thus for large A the RG equation can be 
written in a more conventional form 

which has the solution 

<(A) ~ -\/|- (31) 

(The RG equation can also be integrated exactly, yielding a slightly more complicated 
form.) As A is a dimensionless nonrelativistic energy, this is equivalent to running as 1/p, 
where p is the momentum, a result familiar from EFT. Similarly, a s 2 s runs as A -3 / 2 and 
af as A" 5 / 2 . 
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This procedure is followed through order (r/6) 4 in Ref. ||, which involves a series of 
coupled RG equations for | d , and a^. Beyond lowest order the running 

is no longer automatically scheme independent, but independence can be restored by 
introducing a finite number of higher derivative counterterms in NLO and NNLO. We 
also stress that the analogy between EFT approaches, where the coefficients of contact 
terms are determined by fitting data, and the Talmi integral expansion is very close. In 
NNLO, for example, the three lowest Talmi integrals completely determine the set of 
matrix elements coupling the n = 1 and n = 2 oscillator shells. We can fit our ss bare 
contact potential to those matrix elements, then use the RG equation to run this potential 
to the SM scale, arriving at the same point as our EFT colleagues. It is thus apparent that 
the detailed short-range form of the bare potential is irrelevant to shell-model-inspired 
effective theories: the needed quantities are Talmi integrals over the short-range part of 
the potential. Any change in the detailed radial behavior of the potential that preserves 
the value of the Talmi integrals will also leave the low-energy theory unchanged. 

3.3. Perturbation Theory 

In the above discussion we were careful to carry through the subscript sr on our Talmi 
integrals, denoting the fact that, at the outset, we divide the potential in some smooth 
way into a long-range piece and into a short-range component. It is clearly the short-range 
component - that rather singular part of the bare potential that is unresolved at the scale 
Asm ~ that should be represented by the contact operator expansion. And clearly, while 
we have just solved the short-range effective interactions problem exactly, we are left with 
a gentler long-range potential 

V lr (A 00 ) = V-V sr (A 00 ) (32) 

whose effects have not been evaluated. This potential, if we have defined sr properly, is a 
soft potential that can be handled very well within the SM space: this is where we exploit 
the marvelous technology of the SM. But left over is a small contribution that couples 
the SM and high momentum spaces. How is this treated? 
The procedure is an expansion of the BH propagator 

QH = Q(H + V s ) + Q(V lr - V ) (33) 

treating the second term perturbatively. This corresponds to insertions of the soft poten- 
tial into an infinite series of hard-core scatterings. Thus the lowest order term is 

V?' = V t + V, w ^L— 3 QV, (34) 

the solution of which is just our renormalized contact interaction. The higher-order terms 
could be easily evaluated if we had a convenient expression for the Green's function 
appearing above. Such an expression exists and depends quite simply on the renormalized 
contact interaction coefficients. The expressions can be found in Ref. 0. 

The net result for the deuteron ground state energy is shown in Fig. 2. This calculation 
employs a 10 Tiuj SM space. If one were to try to evaluate the SM effective interaction 
perturbatively, the solid line is obtained. The answer oscillates rather wildly, order by 
order in perturbation theory, failing to converge even after 20 orders of perturbation 
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theory. In contrast, the dashed line is the result of the procedure described here. The 
oscillations damp quickly to an answer than is accurate to a few tens of kilovolts. It 
appears, even in lowest order, that rather accurate results can be obtained perturbatively. 




Figure 2. Calculations of the binding energy of deuterium in a IOKusm included space 
augmented by perturbation theory. The solid result is a standard perturbation treatment 
around the harmonic oscillator potential. The answer is very slowly converging. The 
dashed line is the ET result in which the hard core (defined in this application as the po- 
tential within 0.7 f) is first treated via the LO renormalization group procedure discussed 
in the text: the hard core is summed to all orders analytically, and then the effects of 
the remaining softer potential are evaluated perturbatively. The convergence improves 
remarkably. The corresponding exact result (Aoo = 200) is shown by the dotted line. 



4. OUTLOOK 

The extensions to three- and four-body clusters have not yet been carried out, though 
nothing in the procedures is specific to the deuteron. The RG equations will, of course, 
involve the full set of general three- and four-body translationally invariant contact op- 
erators that can be constructed, just as EFT treatments of the three-body problem were 
required to introduce three-body contact terms. Thus the RG equations become more 
complicated, but not conceptually more difficult. If one began with a bare interaction 
that includes a short-ranged three-body force, the coefficients of the three-body contact 
terms would initially be nonzero, then evolve. If there is no such force, the initial value 
is zero, so that the three-body effective interaction is then entirely induced. 
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If this program can be extended to heavier systems, some remarkable simplifications 
might occur. There exists no LO interaction that involves more than four nucleons: at 
most four nucleons can exist in a relative s-state. Thus, the LO hard-core scattering 
problem in Pb is no more difficult than the corresponding problem in 4 He, in some sense. 
This is an explicit realization of Brueckner's idea of a cluster expansion at high momenta. 
If such a system could be reduced to perturbation theory, as we have achieved in the 
deuteron, some apparently unsolvable problems (such as how to properly handle the 
projection operator Q in such systems) might be overcome. Clearly there are many 
hurdles before one gets to this point, but we are cautiously optimistic that more progress 
will be made. 

In closing, it has been a great pleasure to speak at this meeting is honor of a dis- 
tinguished colleague and friend. We thank Achim Richter and the organizers for the 
opportunity to discuss this work. 
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